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Abstract 

This paper derives for non-linear, time-varying and feedback lineariz- 
able systems simple controller designs to achieve specified state-and time- 
dependent complex convergence rates. This approach can be regarded as 
a general gain-scheduling technique with global exponential stability guar- 
antee. Typical applications include the transonic control of an aircraft with 
strongly Mach or time-dependent eigenvalues or the state-dependent com- 
plex eigenvalue placement of the inverted pendulum. 

As a generalization of the LTI Luenberger observer a dual observer de- 
sign technique is derived for a broad set of non-linear and time- varying sys- 
tems, where so far straightforward observer techniques were not known. The 
resulting observer design is illustrated for non-linear chemical plants, the 
Van-der-Pol oscillator, the discrete logarithmic map series prediction and 
the lighthouse navigation problem. 

These results 123] allow one to shape globally the state- and time-dependent 
convergence behaviour ideally suited to the non-linear or time-varying sys- 
tem. The technique can also be used to provide analytic robustness guaran- 
tees against modelling uncertainties. 

The derivations are based on non-linear contraction theory lITSl . a com- 
paratively recent dynamic system analysis tool whose results will be re- 
viewed and extended. 
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1 Introduction 



Non-Linear contraction theory [Ii[l7l[IS[I3|23|2ni2l|2l|24l[36l[37lisa com- 
paratively recent dynamic analysis and design tool based on an exact differential 
analysis of convergence. Similarly to chaos theory, contraction theory converts a 
non-linear stability problem into a LTV (linear time-varying) first-order stability 
problem by considering the convergence behaviour of neighbouring trajectories. 
Global convergence can be concluded since a chain of neighbouring converging 
trajectories also implies convergence over a finite distance. A brief summary of 
contraction theory is given in section [3l 

Whereas chaos [131 and LTV theory [[Till in section [21 compute numerically 
the transition matrix and hence the time-averaged convergence rates in form of 
the Lyapunov exponents, contraction theory provides explicit analytical bounds 
on the instantaneous convergence or contraction rate. Note that the incremental 
stability approach in Il2| extends these instantaneous analytical contraction rates 
to the integrated convergence of neighbouring trajectories. 

Since contraction theory assesses the convergence of all neighbouring trajec- 
tories to each other, it is a stricter stability condition than Lyapunov convergence, 
which only considers convergence to an equilibrium point. It is this difference 
which enables observer or tracking controller designs, which do not converge to 
an equilibrium point. Also, contraction convergence results are typically expo- 
nential, and thus stronger than those based on most Lyapunov-like methods. 

So far contraction analysis relied in section l3l on finding a suitable metric to 
bound the contraction rate of a system. Depending on the application, the metric 
may be trivial (identity or rescaling of states), or obtained from physics (say, based 
on the inertia matrix in a mechanical system), combination of simpler contracting 
subsystems [181, semi-definite programming [21], sums-of-squares programming 
[[31, or recently contraction analysis of Hamiltoninan systems [(24|. 

This paper fTi\ shows that the computation of the metric may be largely sim- 
plified or indeed avoided altogether by extending the first-order exact differential 
analysis to the placement of state-or time-dependent contraction rats of n-th-order 
(n > 1) continuous systems in controllability form 

x(") = f (x, x("-i), t) + G(x, x("-i), t)u 

with A^-dimensional position x, M-dimensional control input u and time t. In 
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addition a dual observer design is derived for smooth n-th order dynamic systems 
in observability form 

^{n) ^ aS""'\x,t) + af'-')(x,t) + ... + a„(x,i) 

with M-dimensional measurement y(x, t), A^-dimensional position x, A'^-dimensional 
non-linear plant dynamics aj(x, t) and time t. 

A similar method to place state- and time-dependent contraction rates will also 
be derived for the corresponding discrete controllability form 

x'+" = f (x^ x^+"-\ i) + G(x^ x^+"-\ i)u' 

with A^-dimensional position x*, M-dimensional control input u* and time index 
i. In addition a dual observer design is derived for smooth n-th order dynamic 
systems in observability form 

^^+n ^ a;+"-')(x\ z) + a^+"-')(x\ i) + ... + a„(x^ i) 

with M-dimensional measurement y(x*, i), iV-dimensional position x\ A^^-dimensional 
non-linear plant dynamics a^ (x* , i) and time index i. The superscript implies 
now and in the following that the function is mapped j times in the future. 

The following example illustrates the relation of this paper to the design of 
standard LTI controllers and shows that for non-linear, time-varying systems, sta- 
ble convergence is not quantified by the linearized eigenvalues, but by the con- 
traction rates as defined in this paper. 

Example 1.1 : Consider the simplified A/C angle-of-attack dynamics 

a + D{qc, Ma)a + K{qc, Ma)a = G{qc, Ma)u 

with angle-of-attack a, dynamic pressure qc{t), Mach number Ma{t) and control in- 
put u. In a generalization of feedback linearization let us now schedule the complex 
eigenvalues Xi(Ma,qc) and X2{Ma,qc) with Ma{t) and qc{t) to reflect this strong 
non-Unear plant dependence in the A/C controller. This yields the hierarchical or cas- 
caded system with zi= a 

Zl = XlZi + Z2 
Z2 = \2Z2 + Ud{t) 
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which implies the control input u 



a = -D{qc,Ma)a - K{qc,Ma)a + G{qc,Ma)u 
= Ud{t) + (Ai + \2)a - (A1A2 - Ai)a 

where Ai and A2 have to be chosen such that u stays real. The key difference to standard 
gain-scheduling techniques (see e.g. |14|) is the term Ai(Afa, Qc). Only with this term 
exponential convergence with Ai and A2 to the desired trajectory, defined by Ud{t), is 
guaranteed. D 



A major point of this paper will be the extension of example [TTTI to the observer 
and controller design of complex state- and time-dependent contraction rates, con- 
sidering the time-derivatives of the contraction rates to make the analysis correct. 

In section|4] state- and time-dependent contraction rates are "placed", as a gen- 
eralization of standard feedback linearization methods (see e.g. [[TT1l ,|[9l or [31]). 
The generalization is that we can choose state- or time-dependent contraction rates 
Xj(zj,t) to simplify u, to handle only piece- wise controllable systems (under- 
actuacted or intermittently controlled systems), as e.g. in the inverted pendulum 
or in legged locomotion, or simply to achieve state- or time-dependent system per- 
formance. In contrast to standard gain-scheduling techniques (see e.g. [14J, Ii26ll . 
11311 ) global exponential stability guarantees of the state- and time-dependent con- 
traction rates are still given. 

Section \5\ derives a corresponding non-linear observer design. It extends the 
LTV Luenberger observer of constant eigenvalues in [|38ll to higher-order non- 
linear systems with designed state- and time-dependent desired contraction rates. 

The corresponding stability analysis of a given higher-order system is pre- 
sented in section [6l This technique also allows to bound analytically the robust- 
ness of a given controller and observer design with respect to modelling uncer- 
tainties. 

Section |71 [8] and |9] extend the controller and observer design technique to the 
discrete case. We e.g. assess the stability of a non-linear price/demand dynam- 
ics, design an observer for the logistic map problem or derive a simple non-linear 
global observer for the standard bearings-only or lighthouse problem, of navi- 
gating a vehicle using only angular measurements with respect to a fixed point in 
space O. The algorithm is non-linear but very simple. It is new to our knowledge, 
and provides explicit global convergence guarantees. The algorithm's stochastic 



4 



version should serve as a simpler and "exact" alternative to approaches based on 
linearization and the extended Kalman filter, both in the pure bearings-only prob- 
lem and as part of more complex questions such as simultaneous localization and 
mapping (SLAM). 

Concluding remarks are offered in section [TOl 



2 Relation to chaos theory 

Contraction theory and chaos make extensive use of virtual displacements, which 
are differential displacements at fixed time borrowed from mathematical physics 
and optimization theory. Formally, if we view the iV-dimensional position x of 
the system at time t as a smooth function of the initial condition x° and of time, 

X = x(x",t) we get 5x = $(0,t) 5x" with the transition matrix $(0,t) = 
i^(x°(0),0,x(t),t). 

Consider now an A^-dimensional, non-linear, time-varying discrete system 

x*+i = f(x\z) 

The convergence behaviour of neighbouring trajectories is then given by the dis- 
crete virtual dynamics 

5x^+1 = F5x* 

with F = ^(x*, i). The transition of any virtual displacement from to i is then 
given by 

(5x' = $(0,2)(5x° 

with the transition matrix 

$(0,i) = F(-^\..F(-*) (1) 

where the superscript (+j) implies that the function is mapped j times in the 
future. 

Consider now an A^-dimensional, non-linear, time-varying continuous system 

X = f(x,t) 
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The convergence behaviour of neighbouring trajectories is then given by the con- 
tinuous virtual dynamics 

= F(x,t)5x 

with F = ^(x, t). The transition of any virtual displacement from to t is then 
given by 

(5x = <l>(0,t)(5x° 

with the transition matrix 

$(0,t) = I + / F{n)dri + [ F(ri) / ' F{T2)dr2dn + ... (2) 
Jo Jo Jo 

which is equivalent to e-l" ^'^^ for a diagonal Jacobian F. 

The Lyapunov components (see e.g. [|3l) simply correspond to the N'th square 
root of the singularvaluesof$(0,t)or$(0,2). Note that the coordinate invariance 
of this dynamics under smooth coordinate transformations is shown for z, t — )■ oo 
in [|3|. The major problem of chaos theory is that in general the above has to be 
computed numerically. 

What is new in contraction theory is that the transition matrices above can be 
exponentially over/under-bounded in analytical form. This will be shown in the 
following section in Theorem [T] and [2l 



3 First-order contraction theory 

Consider now an A^-dimensional, non-linear, time-varying, complex continuous 
system 

X = f(x,t) 

The convergence behaviour of neighbouring trajectories is then given by the con- 
tinuous virtual dynamics 

OX = 7:— (x, tjox 
ax 

Introducing a general complex A^-dimensional virtual displacement 5z = 0(x, t)5x 
leads to the general virtual dynamics 

d 

—6z = F6z 

dt 



with complex F = ^0 + ©^^© ^. The rate of change of a differential length 
5s = V 5z*'^5z > can now be bounded by 

d _ 5z*^ (F*^ + F) 6z 

Amines < ~JI^^ ~ 7^ — -^max^S 

dt 2os 

where Amax(Amin) is the largest (smallest) eigenvalue of the Hermitian part of F. 

Recall that a complex square matrix A is said to be Hermitian if = A*, 
where ^ denotes matrix transposition and * complex conjugation. The Hermitian 
part of any complex square matrix A is the Hermitian matrix 1/2 (A + A*^) . 
All eigenvalues of a Hermitian matrix are real numbers. A Hermitian matrix A is 
said to be positive definite if all its eigenvalues are strictly positive — this implies 
in turn that for any non-zero real or complex vector x, one has x*^Ax > 0. 

Let us now define a finite distance s = min. f^^._ 5s > between two 
arbitrary trajectories xi and X2 of the dynamics as the minimum path integral 
over all connecting paths s [i25il . The rate of change of a finite length can now be 
bounded by 

AminS < S = min / -i:^^ ^ Xuia.yiS 

Jx{s)=xi 

where Amax(Amin) is the largest (smallest) eigenvalue of the Hermitian part of F 
along the path s. 



The basic theorem of contraction analysis 11181 [1911 can hence be stated as 

Theorem 1 Consider the deterministic system x = f (x, t) , where f is a differ- 
entiable nonlinear complex function of:K within . 

Any trajectory Xi with a distance s = min^ J^^g^^^_^ V 5x*^M(5x > to a given 
other trajectory X2 in a metric M(x, t) = ©(x, t)*^ 0(x, t) > exponentially 
converges to X2 within the bounds 

Amax ^ ~ ^ Amin (3) 
S 

Amax (Aminj is defined as the largest (smallest) eigenvalue of the Hermitian 
part of the generalized Jacobian 

ax / 
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in the ball of radius s around X2. 

The system is said to be contracting (diverging) for uniformly negative Amax (uni- 
formly positive Aminj- The system is said to be semi-contracting (semi-diverging) 
for negative Amax (positive Xmin) (^nd indifferent for Amax = Amin = 0. 

For a u.p.d. and bounded metric also the distance miiio f'A V 5x*^5x con- 

^ * Jx(s)=xi 

verges uniformly exponentially with the rates above, where however initial over- 
shoots can occure. 

Note that the region of convergence of two arbitrary trajecctories with dis- 
tance dynamics s/s in (|2]) can be extended beyond the contracting region with 
Lyapunovs direct method for the specific case that explicite orthonormal Carte- 
sian coordinates z(x, t) with dimension > N exist as 

■ 5Rf(Zi-Z2)*^(Zi-Z2)) 
S (Zi - Z2) (Zi - Z2) 



Note that the theorem above also applies to non-differentiable f if Xmax and 
Xjnin are defined over any limit Um^y^^o-^ instead of the term 

Note that for a semi-contracting system (i.e. with negative semi definite F) 
we can conclude on asymptotic convergence if the indefinite subspace of the sym- 
metric part of F becomes negative-definite in one of the higher time-derivatives 
of Sx^Sx before it eventually becomes positive definite since (5x cannot get stuck 
as long as it is unequal zero. 

It can be shown conversely that the existence of a uniformly positive definite 
metric with respect to which the system is contracting is also a necessary condi- 
tion for global exponential convergence of trajectories. In the linear time-invariant 
case, a system is globally contracting if and only if it is strictly stable, with F sim- 
ply being a normal Jordan form of the system and the coordinate transformation 
to that form. 

The following example shows how for complex systems the contraction region 
of neighbouring trajectories and the region of convergence of trajectories with a 
finite distance can be computed with Theorem [B 

Example 3.1 : Let us now schedule non-linear complex contraction rates for a second- 
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relative radial 




Figure 1: Convergence of two trajectories with finite distance 



order system by requiring the first-order complex dynamics 

z = -^z'^ - 2z + Udit) (5) 

with complex contraction rate A = — z — 2 of Theorem[T] In principle any differentiable 
complex function can be used here to schedule the state-dependent complex contraction 
rates as we want. 

The convergence rate of an arbitrary trajectory zi to another trajectory Z2 is 

S [Zl - Z2)*{ZI - Z2) 2 

according to ^ Theorem [T] This region of convergence is naturally larger then the 
contraction region 3ft(A) < 0. 

The complex dynamics is illustrated in figure [2]for Ud = 0. We can see that the K(A) 
decreases to the right. We find exactly two equilbrium points at zi = and Z2 = —4 
with constant distance - = 0. 

s 

The complex dynamics is with x = and y = ^(z) equivalent to 

X = -^x^ + ^y^ - 2x + u{t) 
y = -xy - 2y 
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Figure 2: Quadratic complex state space dynamics 
with corresponding Jacobian 

/ -X - 2 y \ 
V -V -x-2 J 

that is contracting with K(A) = — x — 2. 

Hence the corresponding real second-order plant dynamics of (|5]) is 

X = -(6 + 3x)x - x^ - 6x^ - 8x + (2x + 4)ud(t) + 
to which the same convergence results apply. D 

For the general A^-dimensional continuous case contraction theory [fTTl [TSll 
can be regarded as time-varying, complex generalization of ifTOl IT3l l30l l33l [T6ll 
with given exponential convergence rate. In addition the introduction of the virtual 
displacements in [fTTl [l8l lead to a generalization of the well-established stability 
and design principles of LTI systems (see e.g. lfT2ll ') to the general non-linear 
and time- varying case. This lead to the practical controller or observer designs in 
|[8l[I7l[l8l[l9l|20l|2ll|23ll2l[33and serves as a basis for this paper. 
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An appropriate metric to show that the system is contracting may be obtained 
from physics, combination of contracting subsystems [18], semi-definite program- 
ming fi2ri . or sums-of-squares programming [4J. The goal of this paper is to show 
that the computation of the metric may be largely simplified or avoided altogether 
by considering the system's higher-order virtual dynamics. 

Similarly, for a discrete system we can state 

Theorem 2 Consider the deterministic system x*+^ = f(x*,i) , where f is a 
smooth non-linear complex function of^^ within . 

Any trajectory x^ with a distance s* = min^i 

to a given other trajectory in a metric M(x*,i) = 0(x', i)*-^ 0(x', i) > 
exponentially converges to Xg within the bounds 

0"max > — > 0"min 0) 

c"max fo'minj IS defined as the largest (smallest) singular value of the general- 
ized Jacobian 

F(x\^) = e(+i)|^e-^ 

in the ball of radius around Xg. 

The system is said to be contracting (diverging) for uniformly negative cTmax — 1 
(uniformly positive (Jmin — 1| The system is said to be contracting (diverging) for 
negative cimax — 1 (positive a^ain — ^) <^nd indifferent for CTmax = o'min = 1- 



For a u.p.d. and bounded metric also the distance minoi ,s . V 5x**^5x* 
converges uniformly exponentially with the rates above, where however initial 
overshoots can occure. 

Note that the region of convergence of two arbitrary trajecctories with dis- 
tance dynamics in ([7|) can be extended beyond the contracting region with 
Lyapunovs direct method for the specific case that explicite orthonormal Carte- 
sian coordinates z\-x\ i) with dimension > N exist as 



s 



_ (zi+^-zn*^(zri-zr^) 



*T (8) 

(zl - zlY^ (zl - zl) 
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This theorem can be regarded as a time- varying, complex generalization of the 
contraction mapping theorem (see e.g. ||6l) to a general metric. This lead to the 
notation Contraction Theory. 

4 Continuous-time controllers 

In this section we consider Vt > a smooth n-th order real dynamic system in 
controllability form 

x(") = f (x, x("-i), t) + G(x, x("-i), t)u 

with A^-dimensional position x, Af -dimensional control input u and time t. The 
controllability conditions under which a general continuous, non-linear, dynamic 
system can be transformed in the form above is well established for feedback 
linearizable systems (see e.g. ^ or ll3T1l ). 

Let us now generalize the well-known LTI eigenvalue-placement in Jordan 
form to the placement of the hierarchical complex dynamics 

Zj = J Aj{zj, t)dzj + 3?(zj+i) for j = 1, ...,p (9) 

with 3? (zi) = X, Zp+i = and where p is given by n minus the number of complex 
contraction rate matrices . Taking the variation of the above implies the time- 
or state-dependent complex contraction rate matrices Aj(zj, t) in 

— 5zj = \j{zj,t)6zj + ^{5zj^i) for j = 1, 

According to Theorem [T] is the stability of this hierarchy given by the definiteness 
of the Hermitian part of Aj . 

Substituting the p dynamics ^ recursively in each other leads to 

Theorem 3 Given the smooth n-th order dynamic system in controllability form 
x(") = f (x, x("-i\ t) + G(x, x("-i), t)u (10) 
with N -dimensional position x, M -dimensional control input u and time t. 
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A controller u that places the complex, integrable contraction rates Aj{zj, t) 
in the characteristic equation 

(|-/V)s^..(|-/A.d)x = (11) 

with 3?(zi) = X and Zj+i = Zj — J Ajdzj implies global contraction behaviour 
with Aj{zj,t) according to TheoremU\ 

p is here given by n minus the number of complex contraction rate matrices 
Aj and 3ft applies to its left-hand term. The open integral j implies a time-varying 
integration constant that can be chosen to shape a desired trajectory in the flow 
field without affecting the contraction behaviour 

The generalization to standard feedback linearization methods (see e.g. BH or 
II3TI ) is that we can choose state- or time-dependent contraction rates Aj(zj , t) to 
simplify u, to handle only piece-wisely controllable systems or simply to achieve 
state- or time-dependent system performance. 

In contrast to well-known gain-scheduling techniques (see e.g. fi^), who also 
intend to achieve state-dependent stability behaviour, we can analytically proof 
global contraction behaviour with Aj{zj,t). Analytic robustness guarantees to 
modelling uncertainties are given in section [6l 

Note that (fTTl) can be modally solved as 

6zj{t) = $(0, t) 3ft(5zj+i(r))<l>(r, 0)rfr + $(0, t)5z° 

with the transition matrix $(0, t) in equation ^ which can be analytically over/under- 
bounded with Theorem[IJ This extends the well-established LTI convolution prin- 
ciple to state- and time-dependent contraction rates. 

Let us first consider LTV systems before we go to the non-linear case: 

Example 4.1 : Consider the second-order real, time- varying dynamics 

X + D(t)x + K(t)x = u(t) 

Real contraction rates Ai{t) and A2(t) imply with the characteristic equation (fTTl) in 
Theorem [3] 

D(t) = -A1-A2 
K(t) = A2A1-A1 
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A complex contraction rate Ai(t) in 



dzi = AiJzi 

implies the real dynamics 

S'k = Re(5x + ImSy 
6y = —Im5x + Re(5y 

with 5x = 5R((5zi), 6y = 9(5zi), Re = K(Ai) and Im = 9(Ai). Rewriting the 
above as second-order dynamics in 5x implies 

D(t) = -2Re - Imlm"^ 

K{t) = -ReRe + Imlm - DRe - Re 

Note that only the additional time-derivative of Ai make this analytic stability result 
correct in comparison to a standard LTI approximation of the LTV system. □ 

Let us now consider real non-linear systems before we go to the complex non- 
linear case: 

Example 4.2: Let us now schedule Xi{zi, Ma, Qc) and X2{z2, Ma, qc) in example 
Example [TTT] (see e.g. ll28l ) in the characteristic equation (fTTI ) in Theorem |3] 



Aofi Q — / Aida = 



d 

m 

with zi = a, Z2 = zi — j Xidzi which is equivalent to 

a = -D{qc,Ma)a-K{qc,Ma)a + G{qc,Ma)u 
= ^ y Xidzi + j X2dz2 

where the time-varying integration constant can be chosen to achieve tracking -behaviour 
of the controller. 

Again the difference to standard gain-scheduling techniques (see e.g. |[T4]| ') is the in- 
tegration over A2 and the time derivative of Ai. Only with these terms exponential 
convergence with the eigenvalues is given. □ 

Let us now go to complex state-dependent contraction rates. This extension 
allows to achieve global stability for partially controllable systems as e.g. the 
inverted pendulum. 
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Example 4.3 : Let us now place for the inverted pendulum without gravity 

X = ncos(x) 

in figure |3] the complex contraction rate A = cos(2;) with z = a; + iy of the complex 
dynamics 

i = sin(z) 

We assume without loss of generality — vr < x < vr. The first-order complex dynamics 
is equivalent to 

X = sin(x) cosh(y) 
ij = cos(a;) sinh(y) 

whose real second-order plant dynamics is 

X = cos(x)x cosh(y) + sin(x) sinh(y)y = cos(x) sin(x)cos/i(2y) 

with the control input u = sin x cosh(2y) that stays bounded for bounded y. 

The chosen contraction rate A is according to Theorem [T] diverging for the lower posi- 
tions cos X > and contracting for the upper positions cos x < 0. 

The convergence rate of an arbitrary trajectory zi to the lower pendulum position Z2 = 
Ois 

s sin(x)x cosh(y) + cos(x) sinh(y)2/ 

~ ~ * ~ 9 \ 9 — ^ (12) 

s zlzi x^ + 

according to Q Theorem [T] We can see that the upper (lower) pendulum position is 
globally stable (unstable) except the trajectory that starts exactly at the lower (upper) 
pendulum position. The corresponding complex dynamics is illustrated in figure |4l 

Let us now choose alternatively the complex dynamics 

z = (1 + i)sin(-) 

with Ai = cos(|). The above is equivalent to 

X = sin(^) cosh(^) — cos(^) sinh(^) 

X y X y 

y = sin(-)cosh(-) +cos(-)sinh(-) 

whose real second-order plant dynamics is 

X = —-cos{x) sinh{y) 

with the control input u = sinh(y) that stays bounded for bounded y. The complex 
dynamics is illustrated in figure [5] We can see that - as designed - every second upper 
position is globally stable / unstable. □ 
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Figure 3: Inverted pendulum 

















3- 


\ 


\l 










: 




\ 


\l 










2- 


\ 


\l 












1- 


\ 


\l 


.\1 












\ 


j/ 






r\ 


















0- 


\ 


1/ 








1 


\\ 




Jv ^ 




















rj 




1 


\\ 






vj 




-1- 


1 


\\ 




P 










1 


\\ 






V 1 


i 




-2- 


1 


\\ 


'l\ 








P 




1 


\\ 


'l\ 


P 








-3- 




'1 











-10 -5 5 10 



X 

Figure 4: Sinus complex state space dynamics of the inverted pendulum 
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Figure 5: Sinus half complex state space dynamics of the inverted pendulum 



5 Continuous-time observers 



In this section we consider \/t >0a smooth n-th order dynamic system in observ- 
ability form 

x(") = aS""')(x,t) + a^"-')(x,t) + ... + a4x,t) 

with M -dimensional measurement y (x, t), A^-dimensional state x, A^-dimensional 
non-linear plant dynamics (x, t) and time t, which is equivalent to 

Xj = Xj+i -I- aj(x, t) for j — 1, n (13) 

with xi = X and x„_,_i = 

Let us now introduce the observer 

kj = Xj+i + aj(x,t) +ej{y,t) - ej{y,t) for j = 1, 

with xi = X and x„+i — that allows to extend the plant dynamics aj with a 
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chosable measurement feedback ej in the equivalent n-th order observer dynamics 

n 

= J2 (a.(x,t) + e,(y,t) - e,(y,t))("-^-) (14) 
i=i 

Let us now generalize the well-known LTI eigenvalue-placement in Jordan 
form to the placement of the hierachical complex dynamics 

Zj = J Aj{zj,t)dzj + 3fJ(zj+i) for j = 1, ...,p — 1 (15) 

with 3?(zi) = X, Zp_|_i = and where p is given by n minus the number of complex 
contraction rate matrices . Taking the variation of the above implies the time- 
or state-dependent complex contraction rate matrices Aj in 

^Szj = \j{zj,t)5zj + 3?(5zj+i) for j = 1, ...,p 

According to Theorem [T] is the stability of this hierachy given by the definiteness 
of the Hermitian part of Aj . 

Substituting the p dynamics (fT5l) recursively in each other leads to 

Theorem 4 Given the smooth n-th order dynamic system in observability form 

= aS"-'Hx,t) + a^"-')(x,t) + ... + a„(x,t) (16) 

with M -dimensional measurement y{^, t), N -dimensional state x, N -dimensional 
non-linear plant dynamics (x, t) and time t. 

An observer 

±j = %+i + (aj (x,t) + ej{y,t) - ej(y,t)) for j = 1, ...,n (17) 

with Xi = X and Xn+i = allows to place with the measurement feedback terms 
Bj the time- or state-dependent, integrable, complex contraction rate matrices 
Aj{zj, t) in the characteristic equation 

= ° (18) 
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with 3f?(zi) = X, 3fJ(zj_|_i) = Zj — f Aj{zj, t)dzj 

The definiteness of the Hermitian part of Aj{zj, t) implies global contraction 
behaviour of the observer state with Aj (zj, t) to the plant state according to The- 
orem\l\ 

p is given by n minus the number of complex contraction rate matrices and 
3? applies to its left-hand term. 

This theorem generalizes the extended LTV Luenberger observer design of 
constant eigenvalues (see e.g. [[TSl. [|26l or [|33|) to non-linear or state-dependent 
contraction rates for non-linear, time- varying systems. It provides a systematic 
observer design technique compared to existing contraction observer designs (see 
e.g. [[D|271[33) 

Note that the global controller in Theorem |3] that uses the state estimates of the 
global observer in Theorem |4] satisfies a separation principle. Indeed, subtracting 
the plant dynamics (fT3] ), eventually extended by a control input G(y, t)u(x, t), 
from the observer dynamics (fTTl ), that is extended by the same control input 
G(y, t)u(x, t), leads with x = x — x and the mid-point theorem to 

• . 9 (a„_j + e„_j) , 

Xj = Xj+i + xforj = l,...,n 

c/x 

with Xri+i = and where ^ is one point between x and x. We can see that the 
Jacobian of the error-dynamics of the observer is unchanged. Since |^ in Theorem 
[3]is bounded the controller represents a hierarchical system [fTSl . As a result is the 
convergence rate of the controller unchanged as well. 

Let us now show how a general Nn dimensional plant 

X = f(x,t) 

with A^-dimensional measurement y = x(x, t) can be transformed to the higher- 
order observability form (fT6l ). A necessary condition is that the mapping 

/ y(x,t) \ 

can be inverted to x(y, y^""^-*, t) such that we get an explicit dynamics (fT6l) 

^(n) =x('^)(x((x,...,x("-l),t),t) 
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Hence a necessary (but not sufficient) observability condition is that the observ- 
ability matrix 

/ \ 
O = : 

with the Lie derivatives ^L°c=^ and U+^c = = Uc^ + j-Vc has 

" " 0x ox ax. at 

piece-wisely full rank. Note that for LTV systems it is also sufficient. 

Let us first consider a linear observer design with time-varying contraction 
rates. 



Example 5.1 : Consider the vertical channel dynamics of a navigation system 

X = a{t) 

with measured altitude y = x and measured vertical acceleration a{t). We want to 
schedule with Theorem|4]the complex eigenvalues Ai(Ma) and X2{Ma) with Ma{t) 
in the observer ([17] ) 



X2 J V '^(*) / V 62 (Ma 

with X = xi to optimize the vertical channel performance for changing altitude mea- 
surement accuracy in sub-, trans- and supersonic. Comparing the second-order ob- 
server error dynamics 

i = (ei(Ma)x) + e2(Ma)x 
at 

with X = X — X to the characteristic equation (fTSl ) 
leads to 

ei(Ma) = A1 + A2 
e2{Ma,Ma) = -X2X1-X2 

The difference to standard gain-scheduling techniques (see e.g. [14]) is the term A2 in 
the feedback gain computation. Only with this term exponential convergence with the 
eigenvalues Ai and A2 is given. 
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Let us consider now observer designs for non-linear systems with time-varying 
contraction rates: 



Example 5.2 : Consider the temperature-dependent reaction A ^ B ina closed tank 

d f ca\ f -I 



dt\ T V -10 



with CA the concentration of A, y = T the measured temperature, and E the specific 
activation energy, where we want to build an observer with designed contraction rates 

Ai(t),A2(t) <0. 

This reaction dynamics is equivalent to the following second-order dynamics in tem- 
perature 

-E -r, E . 

T + ^T^ = -e-TT 
Letting x = e t dT yields the plant in observability form 

X = ai{x) 

E 

with ai{x) = — f e dx. Let us design the observer (fTTl) with estimate x = xi 

\ = ( ^2 + ai \ / ei{y) - ei(y) 
X2 J V J \ e2(y) - e2(y) 

with designed time- varying contraction rates Ai(t), \2{t) < 0. Comparing the equiva- 
lent second-order observer dynamics in x 

x-x = ^ (ai(x) - ai{x) + ei(xi) - ei(x)) + e2(xi) - e2(x) 
with the characteristic equation (fTSl ) 

f4-Ail(x-x)=0 



dt ) \ dt 
leads to the non-linear feedback gains 



ei(y) + ai(y) = (Ai + A2)y 

e2(y) = -(AiA2 + A2)y 



□ 
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The following example gives an explicit equation for the feedback gains of time- 
dependent contraction rates: 

Example 5.3 : Consider the n-dimensional non-linear system dynamics 
= aS"-') (x, t) + 4''^^ (x, t) + . . . + a„(x, t) 

with non-linear plant dynamics aj(x, t) and measurement vector y(x, t). 

Comparing the re-th order dynamics (fT4l) of the observer ([17] ) to the characteristic equa- 
tion dlSl ) of real time- varying contraction rates Aj{t) implies the feedback gains 

n 

ei(y,t) +ai(x,t) = ^A^x 

i=i 

en(y,i) + an(x,t) = adA„-- -{ad A2-Ai)x 
with adoF = -GF - G. □ 

Finally let us consider a non-linear observer with state-dependent contraction 
rates: 

Example 5.4 : Consider the Van-der-Pol oscillator 

X = di{x) + a2(x, t) 

with ai{x) = X — ^ and measured y = x. We want to build an observer (fT/l) with 
estimate x = xi 

\ X2 J V 02(^1, J \e2{y)-e2{y) ) 

with designed contraction rates Ai(a:), A2 < 0. Comparing the equivalent second-order 
observer dynamics in x 

X - X = (ai(x) - ai(x) + eUxi) - ei(x)) + a2(x,t) - a2(x,t) + e2(£i) - e2(x) 
dt 

with the characteristic equation ( fTSl ) 

X — X / (Ai + X2)dx + A2 / Xidx = 

at 
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leads to the non-linear feedback gains 

ei(y) + ai(y) = j\\i + \2)dy 
e2{y,t) + a2{y,t) = -X2 J Xidy 



□ 



6 Continuous higher-order analysis and robustness 

Consider for t > the n-dimensional (n > 1) system 

x(") = f(x,...,x("-i\t) 
with A^-dimensional position x. 

In Theorem [3] and m the characteristic equation of the dynamics above is zero 
since we use the observer or controller feedback to precisely achieve the charac- 
teristic equation. For such a given controller or observer an additional modelling 
uncertainty d may have to be considered on top to the designed characteristic 
dynamics. This introduces the idea of the existence of a perturbation d in the 
characteristic equation if we analyse a given ODE. 

Based on this thought let us approximate this dynamics with the complex, 
integrable contraction rates Aj(zj, t) - that eventually correspond to the designed 
contraction rates - in the distorted characteristic equation 

|-/M)s...(|-/A,<i)x = d(z„.,z„t) 

with 3ft(zi) = X and 3ft(zj+i) = Zj — J Ajdzj. 

Taking the variation of the above we get 

* r 9d ^ dd ^ 

dzp - ApdZp = —dZp + ... + — dzi 

The main idea is to construct in the following an exponential bound on the virtual 
displacement 5zi over p time-derivatives, rather than over the first time-derivative 
as in[18|. 
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Let us first bound the higher-order term part by taking the norm of the above 



|5zp — Ap5zp| < 



dd 



dZr. 



\6zp\ + ...+ 



dd 



dzi 



\Szi\ 



(19) 



where now and in the following the norm of a matrix is the largest singular value 
of that matrix and the norm of a vector is the root of the vector multiplied with its 
conj ungate vector. 



Let us now select a real r](zi{t), ...Zp{t),t) that fulfils 



(20) 



Vt > with V^^r] = Dr] + rfUri > 0, L°r] = 1 and let us bound the initial 
conditions at t = with real and constant K > as 



6zA < K L^-^T] ei'o(''+^— 1< 7 < 



JSP 



(21) 



where Xmaxi'Ziit), ...Zp{t),t) is the largest eigenvalue of the Hermitian part of all 
\j{zj, t). Hence with ([201) and dJl]) we can bound ([l9l) at t = as 



\6ip - \p6zp\ < K e/oC'^+A™^)^* f il^lL^r/ + . . . + ||^|^ 

V 9zp dzi J 



(22) 



Theorem[T]on 5zj = Aj5zj + 'Si{5zj^i) with the bounded distortion (|2T]) and (l22l) 
implies at t = 

\6z,{t + dt)\ < \6z,{t)\ + K {\ma.L^-^v{t) + L^v{t)) e^>'+^-^Ut 
< K U'^r]{t + dt) ei'o^'"(''+^™-), 1 < j < p 

which implies with complete induction that (|2TI) and (l22l) hold Vt > 0. Using the 
above this allow to conclude: 



Theorem 5 Consider for t > the n-dimensional (n>l} system 

x(") = f(x,...,x("-i\t) 
with N-dimensional position x. 
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Let us approximate the above dynamics with the integrable complex contrac- 
tion rates Ajizj, t) in the distorted characteristic equation 

(|-/V)s...(|-/A,.)x = d(z,...z„t) (23) 

with 3fJ(zi) = X and Zj+i = zj — J Ajdzj. 

Bounding the effect of the distortion d with a real i]{zi{t), ...Zp{t),t) that fulfils 
Vt > 

i,,,>||L|i,,+ ...+||L| (24) 

with U^^rj = L^rj+TjUrj > 0, L°rj = 1 andthe largest eigenvalue Xmaxi'^iit), ...Zp{t),t) 
of the Hermitian part of all Aj{zj, t) leads to global contraction behaviour with 

Aj{zj,t)+7] (25) 

according to Theorem \J} where initial overshoots are bounded Wt > with con- 
stant K > by 

\6zj\ < K U-^7]{t = 0) e^o 'J+^-x.rft^ l<j<p (26) 

p is given by n minus the number of complex contraction rate matrices Aj and 
3? applies to its left-hand term. 

Interpreting Aj{zj, t) as the desired contraction rates in Theorem [3] or [Hallows 
to bound the potential instabilities which come from modelling uncertainties of the 
plant. I.e. it allows to prove robustness for modelling uncertainties with a bounded 
de-stabilizing divergence rate rj. Note that additional time- varying errors in the 
control input do not affect the contraction behaviour, but the desired trajectory 

Xd(t). 

If we cannot design Aj(zj, t) per feedback then we have to approximate Aj 
to minimize the distortion d, e.g. by transforming the higher-order system in its 
reduced form (i.e. a form in which f is independent of x*^"^^)). This is illustrated 
in the following examples: 

Example 6.1 : Consider the general (a)periodic dynamics 

.. dil d^U . dV ^ 
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with potentials U{x, t), V{x, t), A^-dimensional position x and where we assume with- 
out los of generality ^^^-^^it) = min ) Vx at a given t. The virtual dynamics 



is 



(5x + 



The distorted characteristic equation (1231 ) is with Ai = — A2 = "^^^lin 



5x + — - + — - Ux = 



^ A 



The remaining instability rj in (1241 ) is then given 

d^U d^U d^V 



r] + rj > 



5x2 ^jj^ Q^2 Q^2 

We can hence bound with Theorem |5]the contraction rates with 



+ 



92[/ d'^U d^V 



9x2 



9x2 ^^2 



Note that for the scalar case with constant damping this condition is equivalent to re- 



quire that the complex poles 
rant of the left-half complex plane. 



92(72 

dx'^ 



92 y 
9x2 



(x, t) lie Vx, t within the ±45" quad- 

□ 



Example 6.2 : Consider the second-order non-linear system 

dV 

x + 2x + — = u{t) 
ox 



with potential energy V = lncosh(x — Xd{t)), that increases the stabilizing force 

dV 
dx 



= — tanh(x — Xd{t)) with the distance to the desired position Xd{t). The corre- 



sponding variational dynamics is 

6x + 26x + 



1 



cosh (x — Xd{t)) 



6x = 



Since the LTI poles lie within ±45 quadrant of the left-half complex plane we can 
conclude with Example |6.1| on contraction behaviour. □ 
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7 Discrete-time controllers 



In this section we consider Vz > a smooth n-th order discrete system in control- 
lability form 

x*+" = f (x\ x^+"-\ i) + G(x\ x^+"-\ i)u' 

with A^-dimensional position x*, M-dimensional control input u' and time index i. 
The controllability conditions under which a general discrete, non-linear, dynamic 
system can be transformed in the form above is well established for feedback 
linearizable systems (see e.g. |fT5ll29ll ). 

Let us now generalize the well-known LTI eigenvalue-placement in Jordan 
form to the placement of the hierarchical complex dynamics 

z}+i = j S,(z},z)dz} + 3ft(z}+i) for J = 1, ...,p (27) 

with 3ft(z 

i) = X*, Zp^i = and where p is given by n minus the number of 
complex contraction rate matrices Sj. Taking the variation of the above implies 
the time- or state-dependent complex contraction rate matrices Sj in 

6z^' = S,(z},z)5z} + 3f?(5z}+i) 

According to Theorem [2] is the stability of this hierachy given by the singular 
values of of S^ . 

Substituting the p dynamics (l27l) recursively in each other leads to: 

Theorem 6 Given the smooth n-th order system in controllability form 

^i+n ^ f (^i^ x^+"-\ i) + G(x\ x^+"-\ i)vi' (28) 

with N -dimensional position x*, M-dimensional control input u* and time index 
i. 

A controller u' that places the complex, integrable contraction rates (zj , i) 
in the characteristic equation 



^) - / Spd ) 3ft... ( (+1) - / Sid ) X* = (29) 
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with = x\3R(z*^J = z*"*"^ — J Sj(z*, z)(iz* implies global contraction 

behaviour with (z* , i) according to Theorem^ 

p is here given by n minus the number of complex contraction rate matrices 
Sj and 3ft applies to its left-hand term. The open integral J implies a time-varying 
integration constant that can be chosen to shape a desired trajectory in the flow 
field without affecting the contraction behaviour 

The generalization to standard feedback linearization methods (see e.g. IfTSl 
|29l ) is that we can choose state- or time-dependent contraction rates Sj (z* , i) to 
simplify u\ to handle only piece-wisely controllable systems or simply to achieve 
a state- or time-dependent system performance. 

In contrast to well-known gain-scheduling techniques (see e.g. [[HI), who also 
intend to achieve state-dependent stability behaviour, we can analytically proof 
global contraction behaviour with Sj(z* , i). Analytic robustness guarantees to 
modelling uncertainties are given in section [9] 

Note that ([291) can be modally solved as 

i 

1=0 

with the transition matrix $(0, z) in equation ^ which can be analytically over/under- 
bounded with Theorem[2l This extends the well-established LTI convolution prin- 
ciple to state- and time-dependent contraction rates. 

Let us first consider LTV systems before we go to real and then complex non- 
linear systems: 

Example 7.1 : Consider the second-order real, time- varying dynamics 

+ D(i)x^+i + K(i)x-' = u{i) 

Real contraction rates 5^i(z) and 5^2(0 imply with the characteristic equation ( [29l ) in 
Theorem [6| 

D(i) = -sS+^)-S2 
K{i) = 
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A complex contraction rate in 

<5z^+i = Si(i)(5zi 

implies the real dynamics 

Jx'+^ = Re^(5x^ + W5y* 
(5y*+^ = -ImMx* + Re*5y* 

with = K(Jzi), 6y' = 9(5zl), Re = 3?(5]i) and Im = Rewriting the 

above as second-order dynamics in Jx' implies 

D(i) = -Re(+^) - Im(+^)ReIm~^ 
K(i) = W+i)ReIm-iRe + Im(+i)lm 

Note that only the change in the time-indices makes this analytic stability result correct 
in comparison to a standard LTI approximation of the LTV system. □ 



Let us now consider the placement of real and state-dependent contraction rates. 



Example 7.2 : Consider the second-order discrete system 
with position and control input u*. 

Let us now schedule ai{z\,i) and a2{z2,'i') with z\ = x^,Z2 = z\^^ — J cridz\ in the 
characteristic equation ( [291 ) of Theorem |6] 

f\^d] [aid]x' = 



This is equivalent to require the control input n* in 

ai{z{+\i + l)dzi+^ + I a2{z\,i)dzl 



where the time- varying integration constant can be chosen to achieve tracking-behaviour 
of the controller to a desired trajectory. □ 



Finally let us consider the placement of complex and state-dependent contraction 
rates. 
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Example 7.3 : Let us now schedule non-linear complex contraction rates for a second- 
order discrete system by requiring the first-order complex dynamics 

/+i = i(z^)2 + nj,(.) (30) 

with complex contraction rate a = \z^\. In principle any differentiable complex func- 
tion can be used here to schedule the state-dependent complex contraction rates as we 
want. 

The convergence rate of an arbitrary trajectory zi to another trajectory Z2 is 

(zl-z|)*^(zl-z|) 2' 1 

according to dH) Theorem |2] This region of convergence is naturally larger then the 
contraction region a <1. 

The complex dynamics is illustrated in figure |6] We can see that \a\ increases from 
the stable origin. We find excactly two equilibrium points at 2:} = and = 2 with 
constant distance = 1. 

The complex dynamics is with = and = equivalent to 



with corresponding Jacobian 



y X 

that is contracting with a = -y/ {x^Y + {y^Y- 

Hence the corresponding real second-order plant dynamics to (l30l ) is 

to which the same convergence results apply. D 



8 Discrete-time observers 



In this section we consider Vi > a smooth ra-th order dynamic system in observ- 
ability form 

^^+n ^ aS+"-^)(x\z) + a(+"-')(x\^) + ... + a„(x\^) 
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Figure 6: Quadratic complex discrete state space dynamics 

with M-dimensional measurement y*(x*, i), A^-dimensional state x*, A^-dimensional 
non-linear plant dynamics aj (x*, i) and time index i, which is equivalent to 

xf 1 = xj^i + a^(x\ for J = 1, n (31) 

with x\ = X* and x^^^ = 0. 

Let us now introduce the observer 

with x.\ = X* and x^_,_i = that allows to extend the plant dynamics with a 
chosable measurement feedback in the equivalent n-th order observer dynamics 

n 

Let us now generalize the well-known LTI eigenvalue-placment in Jordan form 
to the placement of the hieracial complex dynamics 

z}+i = J Aj{z), t)dz) + 3ft(z}+J for J = 1, ...,p - 1 (33) 

31 



with 'Si{z\) = x', Zp^]^ = and where p is given by n minus the number of 
complex contraction rate matrices Sj. Taking the variation of the above implies 
the time-or state-dependent complex contraction rate matrices S* in 

5z}+i = S}(z},2)5z} + 5z}+iforj = l,...,p 

According to Theorem [2] is the stability of this hierachy given by the largest sin- 
gular value of Sj. 



Substituting the p dynamics (1331) recursively in each other leads with to 

Theorem 7 Given the smooth n-th order dynamic system in observability form 

x*+" = ai+"-'^(x\z) + a^+"-')(x\^) + ... + a„(x',^) (34) 

with M -dimensional measurement y(x*,i) of the N -dimensional state x*, A^- 
dimensional non-linear plant dynamics aj(x*, i) and time index i. 

An observer 

x;+^ = x}+i + aj(x\ i) + ej{y\ i) - ej{y\ i)forj = 1, n (35) 

with x\ = X* = 3ft(zi) and x^+i = allows to place with the measurement 
feedback terms ej the time- or state-dependent, integrable, complex contraction 
rate matrices Sj(z*-, i) in the characteristic equation 

(^(+1) - j Epd^ (^(+1) - J Sirf^ x^ = (36) 

with sR(zi) = x\ 3?(z;.+i) = zj+i - / A,(z;.,i)c?z;.. 

The largest singular value o/S* (z* , i) implies global contraction behavior of 
the observer state with (z^ , i) to the plant state according to Theorem^ 

p is given by n minus the number of complex contraction rate matrices E* and 
3? applies to its left-hand term. 

This theorem generalizes the extended LTV Luenberger observer design of 
constant eigenvalues (see e.g. ifTSl . Il26l or [[38l ) to non-linear or state-dependent 
contraction rates for non-linear, time-varying systems. 
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Note that the global controller in Theorem [6] that uses the state estimates of the 
global observer in Theorem |7] satisfies a separation principle. Indeed, subtracting 
the plant dynamics (|3TI) . eventually extended by a control input G(y', i)u*(x*, i), 
from the observer dynamics (l35l) . that is extended by the same control input 
G(y*, i)u*(x*, z), leads with x* = x* — x* and the mid-point theorem to 



(^,^)x'forj = l,...,n 



with X 



n+l 



and where ^ is one point between x* and x*. We can see that 
the Jacobian of the error-dynamics of the observer is unchanged. Since |^ in 
Theorem [3] is bounded the controller represents a hierarchical system ifTSl . As a 
result is the convergence rate of the controller unchanged as well. 

Let us now show how a general Nn dimensional plant 



x*+i = f(x\i) 

with A^-dimensional measurement y* = x*(x*, i) can be transformed to the higher- 
order observability form (|34l ). A necessary condition is that the mapping 



\ 



y^+"-i(x\z) / 



can be inverted to x*(y*, y*"*"" \ i) such that we get an explict dynamics (l34l) 



i+n 



x"+"(x^(y\...,y 



i+n— 1 



Hence a necessary (but not sufficient!) observability condition is that the observ- 
ability matrix 

/ L^c \ 

= 



V 



with the Lie derivatives Bill L°c = |4 and L^+^c = % 



i+j+i 



^i^(+i)|f(+j)has 



piece-wisely full rank. Note that for LTV systems it is also sufficient. 

Let us now consider the observer design of a specific non-linear problem be- 
fore we go to the general non-linear case: 
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Example 8.1 : Consider the logistic map dynamics 

^ ^ 5. (1 _ 

with measured state = U and unknown constant gain c*. We can rewrite the above 
as second-order system 

^,,,2 ^ y+1 y+i (1 - y+i) 
5*(1 - 6^) 

Introducing the complex state x* = In 6* we get 

X — ^1 +0-2 

with ai = ln(6*6*(l - U)) and 02 = - ln(6^(l - 6*)). 
We want to build an observer (1351) with estimate = Xi 



xi+' J { a2(xO J ^ \ e2{f) - e2(,y') 



with designed constant contraction rates |iTi|,|(T2| < 1. Comparing the equivalent 
second-order dynamics (l32l ) in x* 

^i+2_^i+2 ^ _ ^ _ ei(x^))^+'^ + (a2(x*) - a2(x*) + ei(x) - ei 

with the characteristic equation (l36l ) 

(^(+1) - y a2d) ((+^) - y aid) -x')=o 

leads to the non-linear feedback 

e2{f) + a2{y') = -ai02y' 

Note that the observer can be transformed back to the real coordinates h\ = e^i , b\ — 

6 2 as 



ii+l ] ~ L^pe2(y')-e2(y') I 

"2 / \ h\{l-b\) ) 
such that we can compute 6' = \)\ and the estimated unknown gain as 

b\{l-b\) 



□ 
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The following example gives an explicit equation for the feedback gains to achieve 
time-dependent contraction rates: 

Example 8.2 : Consider the n-dimensional non-linear system dynamics 

^^+n ^ aS+"-^)(x\i) + a(+"-') + ... + a„(x\i) 

with non-linear plant dynamics aj(x*, i) and measurement y*(2;*, i). 

Comparing the n-th order dynamics (l32l ) of the observer (1351 ) to the characteristic equa- 
tion (|36l ) of real time- varying contraction rates Xlj(i) implies the feedback gains 

e„(ySi) +a„(x^^) = nj^^S^x^ 

n 

ei(y\i)+ai(x\i) = V 

□ 



9 Discrete higher-order analysis and robustness 

Consider for i > the n-th dimensional (n > 1) system 

with A^-dimensional position x*. 

In Theorem [6] and |7] the characteristic equation of the dynamics above is zero 
since we use the observer or controller feedback to precisely achieve the charac- 
teristic equation. For such a given controller or observer an additional modelling 
uncertainty d may have to be considered on top to the designed characteristic 
dynamics. This introduces the idea of the existence of a perturbation d in the 
characteristic equation if we analyse a given ODE. 

Based on this thought let us approximate this dynamics with the complex, 
integrable contraction rates (z*, i) - that eventually correspond to the designed 
contraction rates - in the distorted characteristic equation 



- 1 Sprf) ((+1) - 1 Sirf) x^ = d(z;, 
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with 9?(zi) = X and 3?(z*_^^) = z*"*"^ — J Sj(z* , i)(izj . For a controller or observer 
of Theorem [6] or |7]d may represent the modelling uncertainities of the system. 

Taking the variation of the above we get 

The main idea is to construct an exponential bound on the virtual displacement 
5x* over p time-steps, rather than over a single time-step as in ifTSl . 



Let us first bound the higher-order term by taking the norm of the above 

- ^pK\ < \d-f\H\ + ■■■ + Ig-rW^-W (37) 

where now and in the following the norm of a matrix is the largest singular value 
of that matrix and the norm of a vector is the root of the vector multiplied with its 
conjungate vector. 

Let us now select a real ri'^{z\, z^, i) > that fulfils 

v'^'-'-V' > l^W'-'-'-V' + ■■■ + \^\ (38) 

CZp oz^ 

Mi > 0. Let us bound the initial conditions at i = with real and constant K > 
as 

\5z]\ < K r^^+^-2...r^^ ^-^7^'= + a^J, 1<J<P (39) 

where cr^naxi'^ii ■■■) z^, i) is the largest singular value of all Sj(Zj, i). Hence with 
(1381) and ^ we can bound dBT]) at i = as 

< Kr^^^^-\..^^W-},U + at,^) (40) 

Theorem [21 on Sz^^ = Aj + ^(zj+i) with the bounded distortion ^ and (gOl) 
implies at i = 1 

|5z;+i| < alJ6z]\+Kv^^^-\..v^ir,-J,{v' + a'^^^) 

< K r^^+^-i...r^*+^ nUo(^' + ^LJ, 1 < J < P 
which implies with complete induction that (|39| ) and (|40|) hold Vi > 0. 

Using the above this allows to conclude: 
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Theorem 8 Consider for i > the n-dimensional (n>l) system 



with N-dimensional position x* at time i. 

Let us approximate the above dynamics with the integrable, complex contrac- 
tion rates Sj(z* , i) in the distorted characteristic equation 

(^(+1) - j Sprf^ (^(+1) - j Sid^ X* = d(z;, ...zi,2) (41) 

with 3f?(zi) = x^ and 3?(z;.+i) = z}+^ - / Sj(z;., i)rfz;.. 

Bounding the effect of the distortion d with a real ri^{z\, ...,Zp,i) > that 
fulfils Vi > 

(9d (9d 

ri'+'^-\..ri' > |T^|r/'+""^.V + • • • + It^I (42) 
azp azi 

and the largest singular value o'lnaxi'^iy ^ip, of all 5]j(z* , leads to global 
contraction behaviour with 

(43) 

according to Theorem |2] where initial overshoots are bounded Vz > with con- 
stant K > by 

\6z]\ < K r^^-\..^°{i = 0) (7?-! + + <ax,), 1 < J < P (44) 

p is given by n minus the number of complex contraction rate matrices Sj and 
3? applies to its left-hand term. 

Figure |7] illustrates the exponential bound (|39| ) which allows short-term over- 
shoots but implies exponential convergence on the long-term. 

Interpreting Sj as the desired contraction rates in Theorem [6] or |7] allows to 
bound the potential instabilities which come from modelling uncertainties of the 
plant. I.e. it allows to prove robustness for modelling uncertainties with a bounded 
de- stabilizing divergence rate if. Note that additional time- varying errors in the 
control input do not affect the contraction behaviour, but the desired trajectory 
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If we cannot design Sj (z* , i) per feedback then we have to approximate 
to minimize the distortion d, e.g. by transforming the higher-order system in its 
reduced form (i.e. a form in which f is independent of x*"*"""^). This is illustrated 
in the following examples: 

Example 9.1 : In economics, consider the price dynamics 

n^+i = f(p\i) 
p^+i = g(n\z) 

with n* the number of sold products at time i and corresponding price p*. 

The first line above defines the customer demand as a reaction to a given price. The 
second Line defines the price, given by the production cost under competition, as a 
reaction to the number of sold items. The dynamics above corresponds to the second- 
order economic growth cycle dynamics 

Contraction behaviour of this economic behaviour with contraction rate 77* can then be 
concluded with equation (l42l ) in Theorem[8]for Si = XI2 = 



5p* 5n* 



(45) 
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That means we get stable (contraction) behaviour if the product of customer demand 
sensitivity to price and production cost sensitivity to number of sold items has singular 
values less than 1. We get unstable (diverging) behaviour for the opposite case. 

Note that this result even holds when no precise model of the sensitivity is known, 
which is usually the case in economic or game situations. Whereas the above is well 
known for LTI economic models we can see that the economic behaviour is unchanged 
for a non-linear, time- varying economic environment. 

The above also corresponds to a game situation (see e.g. ll34l or Q) between two 
players with strategic action and n*. Both players optimize their reaction g and f 
with respect to the opponent's action. We can then again conclude for (1451 ) to global 
contraction behaviour to a unique time-dependent trajectory (in the autonomous case, 
the Nash equilibrium). □ 



Example 9.2 : Consider the general dynamics 



+ 



with potentials [/(x', i), l/(x', i), A^-dimensional position x* and where we assume 



without loss of generality that the singular values of 



_&fu 

<9(x02 



correspond Vx* at a 



given i to the minimal singular values of 



The virtual dynamics is 



5x^+2 + 



5(x 



i\2 



9(x 



i\2 



(5x 



i+1 



5(x 



i\2 



5x* = 



The distorted characteristic equation (|4TI ) is with "Ej 



a(x»)^, 



(+1) 



The remaining instability ry* in (l42l ) is then given by 



9(x02 



<5x' 



5(X^)2^.^5(X02 5(^.)2 

We can hence bound with Theorem |8]the contraction rates with 



d{x 



i\2 



+ 



92 [/ 



Q2U (+1) Q2y 



i\2 



9(x^)2 
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For the scalar case with constant damping this condition is equivalent to require that 
the complex poles ± - ^^{x\i) he \/x\ i within the green square 

of the complex plane in figure [8] D 



0.8- 
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Figure 8: LTI stability circle and non-linear contraction square in complex plane 



Example 9.3: Consider the 2D lighthouse problem in figure |4]of navigating a vehicle 
using only azimuth measurements to a fixed point in space. The dynamic equations 
of the vehicle's motion are 

x*+i = + 

with 2D position x* = {x\,x\)'^ and control input u* = (n^jul)^. The vehicle 
measures only the azimuth to the lighthouse, = tan ip^ = ^ . 

Consider now the observer 

x*+i = x' + u' + (7'' - 1) ^"^sin l' ) ^ ~ ) 

From Theorem |2j this observer is semi-contracting for — 1 < 7* < f. Since the true 
dynamics is a particular solution of the observer dynamics we can then conclude on 
global convergence of x* to x. 
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y, A 



center of vehicle 



angle measurement a , 



fixed point in space 



Figure 9: Lighthouse navigation 

In the case of no model or measurement uncertainty, the optimal choice of 7* is 0. 
Otherwise, the choice of 7* should trade-off the effect of these uncertainties, as e.g. in 
the contraction-based strap-down observer of |[39]| . 

Let us now compute the Jacobian of (l46l) 

5r _ / cos V'* sin \ / 7' \ / cos — sin ^/;* \ 
V — sin V'* cos V'* / V *-* 1 y \ sin V'* cos V'* / 

We can compute for constant 7* = 7 e.g. with MAPLE the square of the largest 
singular value of ^t+t ^ as 

^(7^ + 1) cos(#+i - ^^f + 72 sin(V'^+^ - ^'f 

+^| cos(V''+i - V*)(7' - 1)11 cos(V''+' - ^'firi" - \f + 472)]-^ (47) 

which simplifies for 7 = to | cos(V'*'^^ — V'*)P- Using equation (|42l) in Theorem[8] 
for a\ = a\ = ^ the exponential contraction rate is for 7 = 

r/^+^T/^ > |cos(V'*+^ - V'*)! (48) 

Thus, we can conclude on contraction behaviour over several measurement updates if 
il)^ changes over different i. 

Let us now illustrate the above results with simple simulations in the 2D case with 
position x\ and x\ over the time index i with measurement y*. 
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Figure [TO]shows the motion of a vehicle with constant velocity vector, which is initially 
tangential to the lighthouse. Due to the tangential motion leads the observer (l46l ) with 
7 = to global exponential convergence to the real trajectory with convergence rate 
dH). 

Figure [TT] shows a vehicle with constant velocity vector radial to the lighthouse. The 
observer ( |46l ) with 7=1 achieves global semi-contraction behaviour, i.e. the tangential 
error disappears, whereas the non-observable radial error remains. 

Consider now the 3D lighthouse problem of navigating a vehicle using only azimuth 
V'* and elevation measurements 6** to a fixed point in space. The position dynamics of 
the vehicle's motion is 

with 3D position x* = {x\,x\,x\)'^ and control input u* = {u\,u\,u\)'^ . The vehicle 
measures only azimuth y\ = tanip^ = ^ and elevation 1/2 = tan 6** = 



'■3 



to the lighthouse. 



\/^l^l"t"^2 

The measurement equations can be rewritten in a LTV form in x* as 

xl tan V"* - x\ = 
xl tan C - ^3 = 
x\ tan C - = 

with tanf = ^ and tanC = *'^°^''T^' - 

^ cos ^ cos 

Consider now the observer 



x' + u' + (a'i) 




This dynamics is a superposition of the 2D-lighthouse problem. Hence we can con- 
clude horizontally with (|47] ) or (1481 ) on exponential convergence over several measure- 
ment updates if tp^ changes over different i. The vertical exponential convergence rate 
is then given by the minimum of 1 6* | or | c* | . 



42 



Figure 10: Tangential movement with respect to lighthouse 




Figure 1 1 : Radial movement with respect to lighthouse 
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Since the true dynamics is a particular solution of the observer dynamics we can then 
conclude on global exponential convergence of to x. □ 



10 Concluding Remarks 



This paper derives, for non-linear time- varying systems in controllability form 

x(") = f(x,...x("-i\t) 

simple controller designs in Theorem [3] and [6] to achieve specified exponential 
state-and time-dependent convergence rates. The approach can also be regarded as 
a general gain-scheduling technique with global exponential stability guarantees. 
The resulting design is illustrated for real and complex time- and state-dependent 
contraction rates, inclusive the inverted pendulum. 

A dual observer design technique is also derived for non-linear time-varying 
systems in observability form 

x(") = aS"-'^(x,t)+a("-')(x,t) + ... + a„(x,t) 
x^+" = aS+"-^)(x\z) + a(+"-')(x\^) + ... + a„(x\z) 

, where so far straightforward observer techniques were not known. The resulting 
observer design is illustrated for non-linear chemical plants, the Van-der-Pol os- 
cillator, the discrete logarithmic map series prediction and lighthouse navigation 
problem. 

These results allow one to shape state- and time-dependent global exponential 
convergence rates Aj (zj, t) and Sj (zj, t) ideally suited to the non-linear or time- 
varying system with the generalized characteristic equation 

|-/aw)3^...(|-/Ux ^ 



X* 



with 3?(zi) = X, 3?(zj+i) = Zj - / Aj(zj-, t)rfzj ( 3fJ(zi) = x\ 3?(z;.+i) = z 
/S,(z;,2)rfz;.). 
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Analytic exponential robustness bounds on general non-linear, time-varying 
distortions d on the right-hand side of the characteristic equation are given with 
7] in Theorem [5] or [H Both theorems can also be used to derive analytic state- 
and time-dependent approximated convergence rates for given general non-linear, 
time-varying higher-order systems. 

Note that the general technique of this paper matches the eigenvalue analysis 
for LTI systems. For non-LTI systems additional time derivatives for continuous 
systems and index changes for discrete systems of the contraction rates have to 
be considered. Only with these changes exponential convergence guarantees with 
the time- and state-dependent contraction rates are given. 
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